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ABSTRACT 



The preliminary development of a full-scale Vertical Takeoff and Landing 
(VTOL) Unmanned Air Vehicle (UAV) for the Close-Range mission was 
completed at the Naval Postgraduate School (NPS). The vehicle was based on 
half-scale ducted-fan investigations performed at the UAV Flight Research Lab. 
The resulting design is a fixed-duct, tail-sitter UAV with a canard-configured 
horizontal stabilizer. Major airframe components are used from previous UAVs 
and include the wings from a U. S. Army Aquila and the ducted fan from the 
U.S. Marine Corps AROD. Accomplishments include: 1.) the design and 

fabrication of a carry-through spar and 2.) the design and construction of an 
engine test stand. The carry-through spar was designed using finite element 
analysis and constructed from composite materials. The purpose of the test stand 
is to measure torque, horsepower, and thrust of an entire ducted fan or an 
individual engine. Completion of this thesis will pave the way for future NPS 
research into the growing interest in VTOL UAV technology. 
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I. INTRODUCTION 



A. THE NEED FOR UNMANNED AIR VEHICLES 

The Unmanned Air Vehicle (UAV) began maturity in a warfare role when in 
1960 an aerial camera was mounted to the airframe of a Ryan Q-2C Firebee 
target drone to adapted the vehicle for photographic surveillance missions [Ref. 
1]. At that moment United States military strategists realized that a UAV could 
be placed in an tactical environment to perform high-risk reconnaissance 
missions without endangering more expensive aircraft and aircrew. Over the 
past 30 years UAVs have been adapted into a multitude of roles including 
surveillance, weather monitoring, communications relay, reconnaissance, over- 
the-horizon targeting, and damage assessment. To perform these missions the 
field of UAVs has been developed to cover the entire spectrum from high-speed 
jet vehicles to low-speed hovering configurations. 

A need exists in the Close-Range role to have a vehicle that exhibits Vertical 
Takeoff and Landing (VTOL) capabilities and is also able to rapidly transit to an 
onstation position. The UAV Joint Project Office (JPO) in Washington D.C. has 
contracted technology demonstrators for a vehicle with a takeoff weight of under 
200 lbs., a 50 lb. payload, and a maximum speed of 150 kts. The aircraft must 
be able to takeoff and land within an area 30m x 60m with an obstacle clearance 
of 15 meters. Many variations of helicopter- type UAVs have been developed but 
are limited to horizontal flight speeds of approximately 70 knots. [Ref. 2, p.50] 
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B. THE NATURE OF UAV RESEARCH AT NFS 

The role of the Unmanned Air Vehicle Flight Research Lab (UAV FRL) at 
the Naval Postgraduate School is to lead and support the advances in UAV 
research and technology. The UAV FRL is uniquely able to fill this role by 
providing a wide base of aeronautical knowledge from both the fixed-wing and 
helicopter communities of all branches of the military. Over the recent years the 
UAV FRL has obtained an Aquila airframe and several complete Airborne 
Remotely Operated Device (AROD) vehicles [Ref. 3, p.73] from the Naval 
Command and Control and Ocean Surveillance Center RTD & E Division 
Detachment to serve as a basis for continuing its research. 

From a design standpoint each of these UAVs is very unique. The Aquila 
possessed a strong, light-weight airframe made from advanced composite 
materials; however, it required additional equipment for a catapult launch and 
net recovery. A powerful ducted fan and stability augmentation system enabled 
the AROD to achieve its VTOL capabilities and to sustain hovering flight. 
Reliance on powered lift allowed the AROD to meet its requirement as a 
stationary camera platform but made it unable to achieve transit airspeeds of 
greater than 30 knots. Working with these two concepts, the UAV FRL is 
developing a vehicle to meet the JPO goals by combining the benefits and 
advances from each of these UAV platforms while minimizing the disadvantages 
that each one possessed. This thesis initiated the research on a full-scale proof- 
of-concept vehicle to meet this goal. The research involves the development of a 
new VTOL configuration, the design of a test stand for ground run-up of the 
powerplant, and the design, testing and fabrication of a new wing center 
structure. This full-scale UAV is based on the past research of Ellwood [Ref. 4], 
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Blanchette [Ref. 5], and Brynestad [Ref. 6]. They developed and tested a half- 
scale model to identify problems associated with a vehicle designed to achieve 
both vertical and horizontal flight. The half-scale design encountered difficulties 
in static thrust available, in excess weight, and in the need for a complex tilting- 
duct mechanism. These deficiencies are corrected in the full-scale vehicle by 
using the high-strength, lightweight airframe components from the Aquila and 
AROD airframes and the superior duct design of the AROD to provide efficient 
static thrust. The full-scale model has a new wing center structure designed by 
computer finite element analysis methods and constructed of advanced composite 
materials. 
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II. BACKGROUND 



A. AQUILA 

Development of the U. S. Army Aquila (XMQM-105) began as a program to 
provide a UAV technology demonstrator. The air vehicle was developed under 
Lockheed's proposal, Remotely Piloted Vehicle System Technology 
Demonstrator Program rRPV-STD^ for the U. S. Army . LMSC-D056091, 30 
August 1974. The Aquila was developed for a principal mission of surveillance, 
target acquisition, fire adjustment, and damage assessment in support of Army 
artillery and ground forces. The vehicle is configured as a flying-wing aircraft 
with a wing span of 12 feet, 3 inches. The general arrangement of the Aquila is 
shown in Figure 2. 1 . [Ref. 6] 




Figure 2.1. Aquila three view diagram. 
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Before final production ended, the Aquila had seen many modifications to 
its structural and aerodynamic design. The Aquila airframe obtained by the 
Naval Postgraduate School is from the manufacture period between December 
1974 to December 1977. Table 2.1 lists the characteristics of the Aquila 
configuration [Ref. 7, p.39]. 



Table 2.1. AQUILA CHARACTERISTICS. 



Wing area (projected) 


31.4 feet^ 


Wing span (overall) 


12.35 feet 


Wing area (reference) 


30.3 feet2 


Wing span (reference) 


11.5 feet 


Root chord (WS 13) 


36.73 inches 


Tip chord (WS 69) 


23 inches 


Mean Aerodynamic Chord (MAC) 


32.26 inches 


Aspect ratio (overall) 


4.86 


Taper ratio 


0.58 


Sweep, Leading edge 


28 degrees 


Sweep, .25 chord 


25.2 degrees 


Dihedral, Trailing edge 


4 degrees 


Incidence (WS 13) 


+3 degrees 


Incidence (WS 69) 


0 degrees 


Airfoil section 


modified NACA 23015 


Forward CG location 


20% MAC 


Aft CG location 


22% MAC 


Maximum gross weight 


160 lbs (at sea level) 


Wing loading 


4 lb/ft2 


Power loading 


12 Ib/hp 
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The significant features of the vehicle are: 

• Swept wing. 

• Shrouded pusher propeller. 

• Removable wings for storage. 

• Lightweight Kevlar® construction. 

• Pneumatic rail launch. 

• Vertical barrier net recovery. 

The airframe is constructed of advanced composite materials. Structural 
strength and stiffness are provided through the use of composite sandwich 
construction. This construction technique reduces the need and extra weight 
associated with the conventional method of rib, bulkhead, and stringer structural 
design. The fuselage is designed of varying elliptical cross section which is 
faired very smoothly into the 28° swept wing to provide a very low radar cross 
section. Power is provided with an aft-mounted ducted propeller. The propeller 
shroud provides increased safety to operating personnel during ground launch 
procedures and protection to the propeller during net capture, while also 
providing directional and longitudinal stability. Flight control was accomplished 
through movable elevens on the wing. A major downfall was that the vehicle 
required a large crew of ground personal and a considerable inventory of 
equipment for launch, control, and recovery. 

B. AROD 

The AROD was designed by the Sandia Research Laboratory in 
Albuquerque, New Mexico. This vehicle was designed for a reconnaissance and 
surveillance mission for the United States Marine Corps. It differed significantly 
from the Aquila in that it was capable of vertical takeoff and landing. The 
AROD possessed no flying surfaces and relied solely on powered lift. Control 
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was obtained through four fixed anti-torque vanes and four moveable control 
vanes, all positioned in the propeller wash of the duct. Its main features were; 

• VTOL flight capabilities. 

• Lightweight graphite construction. 

• Compact size. 

• Minimal crew and support equipment. 

• Advanced digital stability augmentation system. 

The AROD suffered from inefficiencies associated with hovering flight for a 
small-diameter powered-lift platform. Since most of the engine output was used 
to maintain powered lift [FIG. 2.2], there was little excess thrust available to 
increase the forward speed. A benefit associated with the AROD was the greater 
static performance provided by the efficient design of the ducted fan. The 
addition of a shroud around the AROD’s three-bladed propeller significantly 
reduced the contraction of the slipstream associated with normal propellers. 
This result increased the mass flow through the fan and produced more static 
thrust than a conventional propeller configuration [Ref. 8, p.l]. Table 2.2 lists 
the measurements taken from the AROD. 

The improvements in thrust for a ducted fan are directly related to its 
geometrical design. Based on the values in Table 2.1, the data from Reference 8 
predict the following benefits in thrust available with the AROD's design. 

1 . A 50% increase in static thrust at a tip speed of 700 fpm over a 

conventional propeller and 66% increase at 800 fpm. 

2. Inlet area ratio of 1.219 is optimal below M=0.3. 

3. Exit area ratio of 1.115 is most beneficial above M=0.1. 

4. The shroud length of 14 inches is optimal below M=0.3. 
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5. The propeller position at 25% chord is the most efficient location for all 

airspeeds. 

6. The three-bladed propeller provides greater thrust over an increased 

number of blades below M=0.1. 

7. The blade tip clearance of 0.00258 (clearance/diameter) optimizes the 

thrust below M=0.2. 



Table 2.2. AROD CHARACTERISTICS. 



Inlet diameter 


29.25 inches 


Propeller diameter, D 


24 inches 


Exit diameter 


26.75 inches 


Inlet area ratio 


1.219 


Exit area ratio 


1.115 


Exterior contour 


tapered rear i 


Propeller location, %c 


25% 


Number of blades 


3 


Tip clearance 


0.03 1 ± 0.005 inches 


Tip clearance (clearance/dia.) 


.00258 


Engine speed, maximum 


8000 rpm 


Engine speed, nominal 


7000 rpm 


Tip speed, maximum 


838 fpm 


Tip speed, nominal 


733 fpm 


Power loading, ^^^ (Pp./p) 


7.25 Hp/ft2 
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Figure 2.2. AROD diagram. 



C. ARCHYTAS TDF 

Initial work on the proof-of-concept VTOL UAV evolved from the theses of 
Blanchette, Ellwood, and Brynestad. This design was named the Archytas after 
the 500 BC Greek scientist and mathematician credited with the design and 
construction of a mechanical flying bird. The Archytas was configured as a 
fixed-wing vehicle constructed from composite materials. It incorporated a 
fuselage and wings that were modeled after a half-scale Aquila. The fuselage had 
the aft center portion removed to accommodate a tilting ducted fan (TDF) 
located at the center of gravity. [Ref. 4] 

The powerplant modeled the AROD by housing an engine-driven fan in a 
circular duct. The fan was manually rotated to direct the thrust downward for 
vertical flight with control given by vanes mounted in the ducted-fan slipstream. 
For horizontal flight the fan could be positioned to redirect the thrust aft. 
Sufficient thrust was available for horizontal flight; however, vertical flight was 
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never achieved. The shroud for the Archytas was designed with an exit area 
ratio of 1.0. Even though this value will prevent contraction of the propeller 
slipstream, it will only increase the static thrust by 26% over a conventional 
propeller [Ref. 8, p.l]. The duct alone was successfully hovered in tethered 
flight with the assistance of a rate gyro [Ref. 6]. 

Three versions of the Archytas were proposed: a tailless model, a short- 
coupled tail design, and a long boom tail configuration [Ref. 4]. The long boom 
Archytas vehicle is shown in Figure 2.3. For purpose of discussion here this 
version will be referred to as the Archytas TDF. Lessons learned from the 
Archytas TDF research were: 

1 . Optimization of ducted-fan design is essential to gain the necessary static 

thrust required for vertical flight. 

2. Incorporation of a tilting mechanism to allow transition between vertical 

and horizontal flight presents a major difficulty. 

3. The structural design of the airframe must be optimized to keep the total 

airframe weight at an absolute minimum. 

The design of the Archytas TDF was strictly intended as a proof-of-concept 
vehicle to examine the flight capabilities in the horizontal and vertical modes. 
The half-scale size did not provide sufficient internal space to house the 
mechanisms necessary for inflight transition, nor was it designed to do so. The 
half-scale platform was conceived to fly in both horizontal and vertical modes, 
but not to perform the transition maneuver. Further development on propeller 
design for increased thrust for the half-scale UAV continues. 
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III. ARCHYTAS TS CONFIGURATION 



A. CURRENT RESEARCH 

This thesis proceeded with the development of the full-scale proof-of-concept 
VTOL UAV. The redesigned Archytas allowed direct assembly of the Aquila 
wings to the AROD duct. The research involved the re-evaluation of the 
Archytas TDF configuration, the construction of a larger and more efficient 
engine test stand, and the design and fabrication of the complex carry-through 
spar required to incorporate the wings and the ducted fan into a single airframe. 
The final Archytas configuration involved many modifications from the initial 
Archytas TDF design. The most significant change included the use of a fixed 
duct rather than a tilting duct. This use enabled a simpler fuselage structure and 
withdrew the requirement for a complex tilting mechanism for the duct. The 
new design incorporated a configuration referred to as a tail-sitter (TS). From 
this concept came the new name Archytas TS. A horizontal stabilizer was added 
to the design to provide longitudinal stability in forward flight. 

B. ARCHYTAS TS STRUCTURAL COMPONENTS 

Configuration of the Archytas TS considered four primary factors. 1.) 
Foremost was the desire that the airframe components from the AROD and the 
Aquila be joined with minimal addition of structural components. This would 
decrease the complexity and lessen the addition of weight. 2.) The new vehicle 
would avoid any modifications to the AROD or Aquila components that would 
reduce their performance or weaken their structural integrity. 3.) The airframe 
design would eliminate any fuselage structure that could interfere with the 
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airflow through the duct. 4.) The landing gear would not add any additional 
weight but still provide stability while the UAV was on the ground. These 
considerations resulted in the design shown in Figure 3.1. 

1. Archytas TS wings 

The wings for the Archytas TS were taken directly from the Aquila 
with no modification to their structure. The airfoil is a modified NACA 23015. 
Coordinates for the airfoil are listed in Appendix A. Aerodynamic loads on the 
wings are supported by two spars positioned at the 25% and 70% chord location. 
The primary material used in the wing construction is a Kevlar/epoxy composite. 
Kevlar® is a high strength, light-weight material with the fiber properties listed 
in Appendix A. [Ref. 7] 

The airfoil shape was maintained without ribs by the use of stiffened 
wing skins. These skins were built using a molded sandwich construction with a 
0.25-inch Nomex® honeycomb core. The wing skins were then bonded to the 
forward and aft spars to form the leading edge d-cell and center torque box. 
The trailing edge of the wing was filled with an epoxy/microsphere slurry to 
reinforce the joint of the upper and lower wing skins. Elevens were located at 
the outboard portion of each wing panel. Construction of the elevens was 
similar to that of the main wings except that the elevens use a solid honeycomb 
core. The elevens pivot on hinges mounted to the aft spar. Actuation of the 
elevens is provided by direct linkage to electronic servos. The laminate 
properties for wing materials are listed in Appendix A. [Ref. 7] 

2. Forward fuselage and canard 

An additional horizontal stabilizer was incorporated in the design to 
provide increased longitudinal stability and control while flying in the horizontal 
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attitude. The canard configuration was chosen to place the horizontal stabilizer 
ahead of the ducted fan. This forward location would keep the tail away from 
the landing ring during vertical takeoff and landing. Configuration of the 
forward fuselage was simplified by using the AROD roll bar. The roll bar 
provides increased protection during takeoff and landing as well as an excellent 
attachment point for the canard. 

3. Landing gear 

Concern for the sizing of the landing gear was given in designing the 
Archytas TS structural configuration. The wing sweep of 28° on the Aquila was 
reduced to 13° to remove the need for extending the landing gear. The 
significant difference in landing ring length is shown in Figure 3.3.1. With the 
smaller wing sweep angle the height of the duct lip was decreased from 51.75 
inches to 35.5 inches. An aerodynamic improvement was also gained with the 
reduced sweep. 

4. Shoulder joints 

To connect the Aquila wings to the AROD duct an additional structure, 
termed shoulder joints for this discussion, was required. These shoulder joints 
were designed to match the end ribs of the Aquila wings and allow for removal 
of the wing sections for transportation and storage. Forward and aft carry- 
through spars were located within the shoulder joints to support forces and 
moments created during flight. The shoulder joints were designed to cover one 
quarter of the duct circumference each to allow two of the four original AROD 
outer duct skins to be re-attached to the upper and lower surfaces of the duct. 
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IV. SPAR DESIGN 



A. DESIGN PROCESS 

A structural modification was necessary to attach the Aquila wings to the 
AROD duct. Four primary factors were considered in the development of the 
spar. 

1. ) Structural support for a 120 lb. airframe at 8 g’s. 

2. ) Minimal addition of any exposed frontal area that would increase profile 

drag. 

3. ) Sufficient faring of the spar to reduce any adverse wing-body 

interference. 

4. ) Composite construction to increase the strength over an equivalent 

weight metal structure. 

The carry-through spar was designed with ring frames [Ref. 9, p.l63] to 
allow the spar structure to pass around the circumference of the AROD without 
interfering with the airflow through the duct. Forward and aft carry-through 
spars were configured to connect with the forward and aft wing spars of the 
Aquila wing. 

The forward spar was sized to fit within the space available between the 
inner and outer skins of the AROD duct. This design completely enclosed the 
upper and lower portions of the spar from any airflow around the duct. The 
inside diameter of the spar was bonded to the inner duct skin to provide greater 
structural integrity. The increased thickness of the duct available at the forward 
location allowed for a thicker spar. The forward spar is 2.5 inches in average 
thickness and has a 1-inch width. The forward spar did not require a larger 
width due to the greater thickness available. The aft spar was designed around 
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the outer extension of the duct ribs that support the control vanes. The AROD 
provided no additional structure to which to bond the aft spar and conceal it 
from the airflow. Therefore, the aft spar was kept to a maximum thickness of 1 
inch from the inner to outer diameter to decrease the frontal area. The spar 
width was increased to 2 inches to account for the smaller thickness. 

Initial design of the spar structure involved constructing a 1/6-scale model 
[Fig. 4.1]. The model allowed visualization of the wing fillet design that would 
be necessary to efficiently blend the wing-duct connection. The typical wing 
fillet with a radius of about 10% of the root chord [Ref. 9, p.l48] was increased 
for the Archytas TS. The fillet in this case was required to house the spar which 
relied on a greater thickness to carry the bending moment around the duct. The 
dimension for the model were selected to maintain the correct wing sweep and 
incidence angle. The model was constructed from a closed-cell blue polystyrene 
foam and the fillet was carved out to create a shape that was modified from the 
Aquila wing airfoil coordinates given in Appendix A. A smooth shape was 
determined that would allow air to travel smoothly from the leading edge to the 
trailing edge without any sharp disruptions in the airflow. The spar positions 
were located on the model and the forward and aft spar cross sections were cut 
out. The spar cross sections were used to determine a planform for the full-size 
spars that could be faired smoothly into the duct. The configuration of the full- 
size spars is shown in Figure 4.2. 

The required dimensions for the full-size spars were calculated and used to 
build a finite element model. The locations of the node points are shown in 
Figure 4.3. Analysis of the spar structure used the MSC/pal2 finite element 
analysis program produced by the MacNeal-Schwender Corporation. The model 
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Figure 4.1. Archytas foam model. 




Figure 4.2. Full size spar configuration. 



20 





Figure 4.3. Forward and aft spar nodal point locations. 
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file and the nodal point locations are given in Appendix B. MSC/pal2 contains a 
composites module for analyzing quadrilateral elements constructed from 
composite laminates. Classical lamination theory is used by MSC/pal2 to 
compute the results for a laminate with up to 200 plies and ten different 
materials [Ref. 10, p. 1.1]. 

B. MATERIAL SELECTION 

The spars were built using moldless composite sandwich construction. This 
technique used a foam core with a composite laminate hand laid on the shear web 
and spar caps. 

1. Core 

The spar core was made from 2 Ib./cu. ft. urethane foam. Urethane 
was selected because it is relatively inexpensive, extremely easy to cut and shape 
and bonds well with epoxy matrix adhesive. The foam construction considerably 
minimized the waste generated during the core fabrication. The cores for both 
full length spars and the half-span test spars were built from segments cut from a 
single 4’ x 8' sheet of 1 -inch- thick urethane foam [Fig. 4.1.1]. This involved 
cutting single pieces for the forward spar and two pieces for each position on the 
aft spar to built up the 2-inch width. Additional foam blocks were cut to 
construct a half span forward and aft spar for load testing. The foam segments 
were cut from the urethane sheet with a hacksaw blade. The edges were planed 
square to produce a good glue joint. The foam pieces were glued together with a 
micro-slurry mixture of 50% microspheres and 50% epoxy resin [Ref. 11, p.36]. 
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Figure 4.1.1. Foam core cutouts. 



2. Adhesive 

The adhesive used to bond the composite fibers together and attach them 
to the core is technically referred to as the matrix. The matrix used on the 
Archytas TS spar is Safe-T-Poxy® two-part epoxy manufactured by Hexcel 
Corporation. A ratio by weight of 100 parts resin to 48 parts hardener was used 
to mix the epoxy. Cure time for the epoxy is 10 hours to touch and 24 hours 
before sanding. The material properties for epoxy are listed in Table 4.2.1. 

Table 4.2.1 MATRIX PROPERTffiS. 



Matrix 


Ei,psi 


E 2 , psi 


'On 


Tensile strength 


% elongation 


Epoxy 


0.62 X 106 


0.62 X 106 


0.34 


8 X lO^psi 


6.0 



3. Reinforcement fiber 

Three types of fiber were used to construct the laminates: E-glass, S- 
glass and graphite. E-glass is an older form of fiber that is commonly found is 
woven fiberglass cloths. The shear webs were constructed from an Owens- 
Coming woven E-glass cloth. The fabric is a 7725 twill weave that is commonly 
used in homebuilt aircraft. Twill was chosen because it has greater strength 
properties over a plain weave for the same fabric weight and amount of fibers. 
This is because material properties are lost by curving of the fibers in the 
weaving process. Figure 4.2.1 compares a twill weave with a plain weave fabric. 
It can be seen that fibers woven in a twill do not cross over each other as often as 
in a plain weave for a given length of fabric. 

The unidirectional fiberglass used in the Archytas TS is Orcoweb S-500 
manufactured by Orcon Corporation. It serves as the primary structural fiber in 
the spar caps because its strength properties are all oriented in a single direction. 
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Unidirectional fibers do not have any weaving crimps and have not been 
weakened by abrasive weaving processes. This material is constructed with 
Owens-Coming S-2® glass. Orcoweb S-500 fibers are held together by a single 
cross fiber at 1.5-inch spacing that is bonded with a thin layer of adhesive. S- 
glass is stronger and stiffer and weighs less than E-glass and is ideally suited for 
low-cost, high-strength, aerospace applications. [Ref. 13] 
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Figure 4.3.1. Types of weaves. 

Graphite was used at the locations on the spar that required greater 
material properties than those available with S-glass. Graphite provides 
maximum stiffness (high modulus) and extremely high compressive and tensile 
strengths at a very low weight. The use of graphite falls into area of advance 
composites. The graphite used on the Archytas TS is manufactured by Hexcel. 
The fibers are held together by a widely space cross woven glass fiber. The 
material properties of E-glass, S-glass, and graphite are shown in Table 4.3.1. 
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Table 4.3.1 FIBER PROPERTIES. 



Fiber 


El, X 106 


E2 X 106 


Gi2 X 106 


1)12 


Tensile strength 


% elongation 


E-glass 


10.5 psi 


10.5 psi 


4.37 psi 


0.2 


250 X lO^psi 


2.4 


S-glass 


12.6 psi 


12.6 psi 


5.17 psi 


0.2 


360 X lO^psi 


2.9 


Graphite 


32 psi 


2 psi 


2 psi 


0.2 


360 X 10^ psi 


1.1 



4. Laminate design 

The laminates were designed using classical lamination theory [Ref. 12, 
pp.l47]. Laminate requirements were identified from the loads to be applied to 
the spar during flight at 8 g's with a minimum safety factor of two. Forces and 
moments were related to strains and curvatures by 
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where N and M are the loads and moments, A, B and D are the constituent 
material stiffness matrices, and e and k are the strains and curvatures [Ref. 10, 
p.6-6]. The A, B and D values for a laminate constitute what is referred to as the 

'ABD matrix'. The values for the ABD matrix were obtained from the 
transformed reduced stiffness matrix, for each laminate. 
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n + 1 

A, = I(Q.).(K-K_.) 

k=l 

1 

B«=Tl(Q,i),(hJ-hL) 

^ k=l 

1 

d„= 5 E(QsU»;- 1 ' 1 ,) 

^ k=l 

The transformed reduced stiffness matrix took into account the ply angle for 
each laminate [Ref. 10, p.6-4] where 

Qi 1 = Qi 1 cos^ p + 2(Q,2 + 2 Qs6 ) sin^ P cos^ P + Q22 sin^ p 
Q22 = Qii sin^ P + 2 (Qi2 + 2 Qg 6 )sin^ Pcos^ P + Q22 cos^ p 

Qi 2 = (Qn + Q22 - 4 Q 66 )sin^ Pcos^ P + Qi2(sin^ P + cos^ p) 

Q 66 = (Qii + Q22 - 2 Qi2 - 2 Qg 6 )sin^ Pcos^ p + Q«(sin^ p + cos^ p) 

Qi6 = (Qi 1 - Qi 2 - 2Q66 ) sin p cos^ p + (Q 12 - Q 22 + P cos p 
Q26 = (Qii - Qi2 - 2 Qfi 6 )sin^ p cos P + (Q,2 - Q22 + ^Q 66 ) sin P cos^ p 



The reduced stiffness matrix, Qij, for each laminate was derived from 
the ply material properties [Ref. 10, p.6-3]. 



Qn = 
Q22 ~ 



El 



(1-0,2=^ 021) 
£2 

(1--U,2*-U2i) 

(l-^:2*^2:) 

Q6. = G 
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The ply material properties are calculated from the fiber and matrix properties 
that create the laminate [Ref. 10, p.6-13]. 
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V Ein, (l Fv) 






E /E - 1 ^ 

l+ ^2f/Ji2n ^ i ,p 

E,,/E,, + l 



E2f/E2. + 1 

^12 “ ^IfEv + ’^lm(^ ~ Ev) 



Gi2 = 



_ ^12ni(^12f Ey) + G]2m(l Ey)) 



0,„(1-Fv) + G,j„(l-Fv) 



where Fy is the fiber volume of the laminate and the subscripts m and f refer to 
the matrix and fiber properties, respectively. 

5. Shear web 

Laminate design was similar for the shear web for both the forward and 
aft spars. Any coupling between extension and bending of the shear web was 
not desired. If coupling occurred, any shear force applied to the spar would 
result in a twisting of the shear web laminate. To avoid coupling the shear web 
was formed from layers that were symmetric about the laminate midplane with 
opposite signs of ply angle orientation, ±(3. This simplified the ABD matrix by 
setting Ai 6=A26=0 and [Bij]=0 [Ref. 12, p.l65]. 
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To minimize the laminate shear strain, £xy , 


it was desired to maximize the 



A^6 value in the ABD matrix. This calculation required the determination of the 
optimum ply angle, |3, which would result in the largest values of in the 

transformed reduced stiffness matrix. With Qn, Q 22 , Q 12 and Q66 values fixed 
from the material properties, the angle, (3, was varied from -90° to -i-90°. Figure 
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4.5.1 reveals that the maximum values were obtained for (5 = ±45°. The program 
used to determine the ply angle optimization was written in MATLAB®. The 
code is given in Appendix C. This ply angle agrees with the industry standard to 
orient material properties at ±45° to the longitudinal axis to provide the highest 
shear properties [Ref. 14, p.386]. The 7725 twill weave fiberglass was ideal for 
application to the shear web. The fabric oriented at ±45° created a regular 
symmetric angle-ply laminate that had optimal shear properties. An individual 
thickness, t, of woven fabric was modeled in MSC/pal2 by placing four layers of 
unidirectional material of thickness t/4 at ±45° orientation. The ply properties, 
laminate layup, and laminate properties are listed in Appendix C. The laminate 
thickness for the shear web varied at different locations across the spar. The 
thickness was chosen to minimize the stress concentration caused from the flight 
loads. The laminate layup for the shear web is shown in Table 4.5.1. The 
location of each laminate corresponds with the placement shown in Figure 4.5.2. 




Ply angle, degrees 

Figure 4.5.1. Shear web optimization. 
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Table 4.5.1. SHEAR WEB LAMINATE THICKNESS. 



Spar 


Location 


Layers 


Forward 


A 


3 


B 


4 


C 


5 


Aft 


D 


4 
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Fwd. spar 

Shear web laminate location Spar cap laminate location 

(4lyp.) A aA (2typ.) 




Aft spar 

Shear web laminate location Spar cap laminate location 



(4 typ.) 



(2 typ.) 




Figure 4.5.2. Spar laminate location. 
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6 . Spar cap 

The laminates for the spar caps were similar in design, but differed in 
material composition. The forward spar caps were required to carry more of 
the flight loads. Additional stiffness and strength were designed into the forward 
spar caps by increasing the laminate thickness and adding more layers of carbon 
fiber at critical locations. Bending stiffness is controlled by the Dy values in the 
ABD matrix. The carbon fibers were located at the outer layers of the laminate. 
This location provided the greatest effect of (h^ -hk_i) on creating larger values 

of Djj when looking at the entire spar height as a single laminate. 

The optimum ply angle was determined in a similar manner as that for 
the shear web. For the spar caps the flight loads would create an Nx value in the 
ABD matrix. This force caused compression on the upper spar cap and tension 
on the lower spar cap. A nonsymmetric laminate with multiple specially 
orthotropic layers was created by placing the graphite on the outside plys with all 
of the principal material directions aligned with the laminate axis [Ref. 12, 
p.l62]. The individual spar caps were not considered symmetric because of the 
difference in material properties for graphite and S-glass; however, the upper 
and lower spar caps with the foam core in between created a symmetric spar. 
The ABD matrix for an individual spar cap took the form of 
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To minimize the normal strain, , due to , large values of An and A\j were 
desired. Figure 4.6.1 shows that the maximum values were obtained for P=0°. 

The laminate thickness for the spar caps varied at different locations 
across the spar. The thickness was chosen to minimize the stress concentration 
caused from the flight loads. The laminate layup for the spar caps is given in 
Table 4.6. 1 The location of each laminate corresponds with the placement shown 
in Figure 4.5.2. Unlike the shear web, the spar cap layers were modeled directly 
in MSC/pal2 as they appear in the actual laminate. 




Ply angle, degrees 

Figure 4.6.1 . Spar cap ply angle optimization. 



Table 4.6.1. SPAR CAP LAMINATE THICKNESS. 



Spar 


Location 


S-glass layers 


Graphite layers 


Forward 


AA 


4 


2 


BB 


8 


3 


CC 


3 


0 


Aft 


DD 


3 


1 


EE 


4 


2 
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V. TEST STAND 



A. TEST STAND REQUIREMENTS 

A fundamental requirement in the calculation of aircraft performance is 
accurate knowledge of powerplant characteristics. For a piston-driven engine, 
accurate knowledge of available horsepower is essential in calculating horizontal 
flight performance. VTOL-configured aircraft are also concerned with vertical 
flight performance, similar to helicopters, and require information on maximum 
thrust available. Procedures used in the power plant analysis for the Archytas 
TDF are described in Reference 5. The test stand for the TDF design consisted 
of two main components. One portion allowed for determination of the 
horsepower by measuring the torque created by the engine over a range of RPM. 
Static thrust output for the engine was obtained by a separate part of the test 
stand by attaching the entire ducted fan to a sliding tray. The sliding tray pulled 
on a spring scale which measured static thrust directly. In principle the 
equipment was sufficient; however, the following deficiencies were observed: 

1. Inability to measure thrust and horsepower simultaneously. 

2. Airflow blockage resulting from the thrust stand configuration. 

3. Friction losses due to poor bearing surfaces for the sliding tray. 

4. Inadequate size in both the torque and thrust components to support the 

larger powerplant for the Archytas TS. 

A recommendation resulting from the Archytas TDF research was to design 
and construct a larger engine test stand that could accurately measure and record 
thrust, torque, and engine RPM [Ref. 6, p.64]. 
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B. TEST STAND DESIGN 

Initial considerations for a larger test stand began simply with a scaled-up 
version of the smaller configuration used for the TDF design. However, the 
mechanical complexity involved in constructing an efficient sliding tray to 
support an apparatus for measuring both torque and thrust was not desirable. 
The final focus was directed at a linear bearing arrangement. A linear bearing 
consists of numerous rows of ball bearings mounted axially along the inner 
circumference of the outer bearing race. These rows of ball bearings allow a 
highly polished, hardened steel shaft to traverse linearly through the bearing and 
rotate simultaneously. An arrangement of two linear bearings mounted in series 
along with a 16-inch shaft was salvaged from a previous NFS experiment. The 
equipment was restored and modified to create a test stand on which to mount the 
entire AROD duct. The test stand was mounted to a large base to provide 
support and stability during use. Wheels mounted underneath the base could be 
extended to roll the entire assembly; for engine tests, the wheels could be 
retracted to place the frame firmly on the ground and prevent any movement. 
The complete configuration of the test stand is shown in Figure 5.1. 

C. TEST STAND ANALYSIS 

Static analysis was performed on the load-carrying components of the test 
stand listed in Table 5.1 to ensure that they could withstand the loads applied 
with the duct in place and additional forces produced during engine run-up. 

The analysis was performed using a strength-of-materials approach [Ref. 15]. 
The maximum normal and shear stress was calculated for each component. 
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Table 5.1 COMPONENT ANALYSIS. 



Component 


Analysis 


Bearing 


Static load. 


Shaft 


Bending, shear, torsion. 


Extension 


Bending, shear, torsion. 


Bracket 


Bending, shear, torsion. 




Figure 5.1. Engine test stand. 
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The material yield stress and the maximum stress values were applied to Von 
Mises' yield criterion to ensure that no component of the test stand would reach a 
point of yield. Calculations for the analysis of each component is given in 
Appendix D. The material yield stress was compared to the Von Mises' yield 
criterion to produce a factor of safety. 

_ - - Material yield stress 

Factor of safety = 

Von Mises' yield stress 

The results from the analysis in Appendix D are listed in Table 5. 



Table 5.2 FACTOR OF SAFETY. 



Component 


Material Yield Stress 


Von Mises' Yield Stress 


Factor of Safety 


Bearings 


1560 lbs.* 


314.5* 


5.0 


Shaft 


275,000 psi 


6358.1 psi 


43.3 


Extension 


36,000 psi 


3479.3 psi 


10.3 


Bracket 


37,000 psi 


1008.2 psi 


36.7 



Bearing normal and actual load ratings. 



D. ENGINE TEST RESULTS 

Engine tests were performed with the AROD mounted on the test stand. A 
6-inch moment arm [App. E] was clamped to the shaft to regulate the rotation. 
Torque was measured with a spring mounted on the end of the moment arm. 
The thrust from the engine pushed on the shaft and exerted a load on an spring 
load cell. The throttle setting was controlled with an electric servo. Engine 
RPM was set with a computer controlling the throttle servo. Torque and thrust 
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readings were obtained simultaneously over the range of RPM listed in Table 
5.3. Horsepower was calculated from the measured RPM using 



Hp = Torque (in.-ibs.) •{2n RPM) • 



1 u ^ 

1 Hp 

550 ~^^ 



\60tamJ \12 ft. / 



see. / 



The test data are plotted in Figures 5.2 and 5.3. 



Table 5.3. ENGINE TEST DATA. 



Computer 


Engine RPM 


Thrust 


Moment 


Horsepower 


Throttle setting 




Ibf. 


Ibf 




50.0 


1929.0 


8.0 


-4.0 


-0.735 


55.0 


2115.0 


10.0 


-3.0 


-0.604 


60.0 


2210.0 


13.0 


-3.0 


-0.631 


65.0 


2850.0 


19.0 


-1.0 


-0.271 


70.0 


3720.0 


25.0 


0.0 


0.000 


75.0 


4050.0 


32.0 


2.0 


0.771 


80.0 


5160.0 


56.0 


6.0 


2.947 


85.0 


5820.0 


62.0 


10.0 


5.540 


90.0 


6030.0 


65.0 


8.0 


4.592 


95.0 


6300.0 


70.0 


8.0 


4.798 


100.0 


6630.0 


73.0 


6.0 


3.787 


105.0 


6960.0 


78.0 


9.0 


5.963 


110.0 


7080.0 


79.0 


11.0 


7.414 


115.0 


7110.0 


83.0 


16.0 


10.830 


120.0 


7380.0 


88.0 


20.0 


14.051 


125.0 


7380.0 


90.0 


18.0 


12.646 


130.0 


7350.0 


90.0 


14.0 


9.796 
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Thrust Reading, T,lb Horsepower 




RPM 

Figure 5.2. Horsepower vs. RPM. 




RPM 

Figure 5.3. Thrust vs. RPM. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

A full-scale VTOL Unmanned Air Vehicle was developed based on previous 
research on a half-scale UAV. The configuration was modified considerably 
from the half-scale design. The airframe was designed to incorporate the wings 
from a U.S. Army Aquila and the ducted fan from a U.S. Marine Corps AROD. 

The goal of this research was to select a new configuration that did not 
require a complex duct-tilting mechanism to transition from vertical to 
horizontal flight. The investigation resulted in a fixed duct, tailsitter vehicle 
with a canard configured horizontal stabilizer. 

Forward and aft carry-through spars were developed to connect the wings to 
the ducted fan. The spars were designed using computer finite element analysis. 
The material selected for the spar construction consisted of a urethane foam core 
and fiberglass/epoxy reinforced laminated. Carbon fiber was added to the 
composite in areas that require increased strength and stiffness. Laminate design 
was selected using finite element analysis. 

An accurate knowledge of engine parameters was required to effectively 
begin development of the UAV. A large-scale engine test stand was designed 
and constructed to obtain powerplant thrust and horsepower measurements. The 
test stand provided attachment for the entire AROD and was supported by a 
linear-bearing arrangement. Engine tests that were performed indicated that a 
more accurate load cell for the thrust and torque measurements was required. In 
addition, the fixed anti-torque vanes interfered with the horsepower 
measurements. 
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B. RECOMMENDATIONS 

The AROD roll bar should provide an ideal attachment point for the canard. 
The canard sizing should be calculated based on the fixed length of the roll bar. 

The strength of the composite layup is very dependent on the construction 
procedures and laminate orientation, The manufacturing process can introduce 
excess voids or resin-rich areas that will lower the composite strength from the 
design limit. Load testing of the half-span test spar is vital to the design of the 
full-size carry-through spar to insure that the finite element design was sufficient 
to size the spar. 

The negative horsepower measurements identified the influence that the duct 
had on the engine torque measurements. To obtain accurate horsepower 
readings it will be necessary to construct an attachment to mount just an engine 
on the test stand. A accurate strain-gauge load cell should be acquired to obtain 
accurate test stand measurements. 



41 



APPENDIX A 



Table Al. WING AIRFOIL COORDINATES*. 



Station 


Upper Surface 


Lower Surface 


0.0 


— 


0.0 


1.25 


3.34 


1.54 


2.5 


4.44 


2.25 


5.0 


5.89 


3.04 


7.5 


6.90 


3.61 


10 


7.64 


4.09 


15 


8.52 


4.84 


20 


8.92 


5.41 


25 


9.08 


5.78 


30 


9.05 


5.96 


40 


8.59 


5.92 


50 


7.74 


5.50 


60 


6.61 


4.81 


70 


5.25 


3.79 (modified) 


80 


3.73 


2.43 (modified) 


90 


2.26 (modified) 


1.00 (modified) 


95 


1 .69 (modified) 


.0239 (modified) 


100 


1.00 (modified) 


.016 ± 0.5 (modified) 



* [Ref. 6, p.82] 
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Table A2. KEVLAR® MATERIAL PROPERTIES.* 





Tension 


Compression 


Material 


Strength 

(ksi) 


Modulus 

(psi) 


Strength 

(ksi) 


Modulus 

(psi) 


Kevlar 

Unidirectional 


170 


10.1 


40 


10.1 


Style Fabric 120 


60 


4.4 


23 


3.5 


Style Fabric 281 


70 


4.4 


26 


4.0 



* [Ref. 6, p.94] 



Table A3. WING LAMINATE PROPERTIES* 



Member 


Material 


Load 


ctmax (psi) 


c^ult (psi) 


Factor of Safety 


Wing 

skins 


Honeycomb 

sandwich 


Compression 


8880 


12,000 


1.35 


Elevens 


Honeycomb 

sandwich 


Compression 


2534 


12,000 


4.74 


Fwd. 

spar 


Kevlar 


Compression 


17,868 


25,700 


1.443 


Aft Spar 


Kevlar 


Compression 


14,758 


25,700 


1.763 



* [Ref. 6, p.69] 
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APPENDIX B 



A. MSC/PAL2 MODEL FILE 

TITLE ARCHYTAS SPAR(MIN) 

NODAL POINT LOCATIONS 3 

1 12.25 -90 0.25 THROUGH 73 12.25 90 0.25 STEP 6 

2 12.5 -90 0.75 THROUGH 74 12.5 90 0.75 STEP 6 

5 12.59375 -90 1.75 THROUGH 77 12.59375 90 1.75 STEP 6 

6 12.65625 -90 2.25 THROUGH 78 12.65625 90 2.25 STEP 6 

3 15 -90 0.75 THROUGH 21 15 -45 0.75 STEP 6 

57 15 45 0.75 THROUGH 75 15 90 0.75 STEP 6 

4 14.96875 -90 1.75 THROUGH 22 14.96875 -45 1.75 STEP 6 

58 14.96875 45 1.75 THROUGH 76 14.96875 90 1.75 STEP 6 

101 14.875 -90 18.75 THROUGH 149 14.875 90 18.75 STEP 4 

102 14.875 -90 20.75 THROUGH 150 14.875 90 20.75 STEP 4 

103 15.875 -90 20.75 THROUGH 115 15.875 -45 20.75 STEP 4 

139 15.875 45 0.75 THROUGH 151 15.875 90 20.75 STEP 4 

104 15.875 -90 18.75 THROUGH 116 15.875 -45 18.75 STEP 4 

140 15.875 45 1 8.75 THROUGH 152 15.875 90 18.75 STEP 4 
NODAL POINT LOCATIONS 1 

33 17 -5.1875 0.75 THROUGH 45 17 5.1875 0.75 STEP 6 

34 17 -5.1875 1.75 THROUGH 46 17 5.1875 1.75 STEP 6 

79 20.25 -3.4583 0.75 THROUGH 83 20.25 3.4583 0.75 STEP 2 

80 20.25 -3.4583 1.75 THROUGH 84 20.25 3.4583 1.75 STEP 2 

85 23.75 -2.25 0.75 THROUGH 89 23.75 2.25 0.75 STEP 2 

86 24.0 -2.25 1.75 THROUGH 90 24.0 2.25 1.75 STEP 2 

27 13.5833 -8.25 0.75 THROUGH 51 13.5833 8.25 0.75 STEP 24 

28 13.5833 -8.25 1.75 THROUGH 52 13.5833 8.25 1.75 STEP 24 

123 16.875 -4.1875 20.75 THROUGH 131 16.875 4.1875 20.75 STEP 4 

124 16.875 -4.1875 18.75 THROUGH 132 16.875 4.1875 18.75 STEP 4 

153 22.75 -1.626 20.75 THROUGH 157 22.75 1.625 20.75 STEP 2 

154 22.75 -1.625 18.75 THROUGH 158 22.75 1.625 18.75 STEP 2 

159 28 -1.1875 20.75 THROUGH 163 28 1.1875 20.75 STEP 2 

160 28 -1.1875 18.75 THROUGH 164 28 1.1875 18.75 STEP 2 

119 13.9375 -8 20.75 THROUGH 135 13.9375 8 20.75 STEP 16 

120 13.9375 -8 18.75 THROUGH 136 13.9375 8 18.75 STEP 16 
301 21.75 .5 -7 THROUGH 302 21.75 -.5 -7 

307 30.428 .5 28 THROUGH 308 30.428 -.5 28 
309 47.482 .5 0 THROUGH 310 47.482 -.5 0 
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311 47.482 2.125 6.263 THROUGH 312 47.482 -2.125 6.263 
313 47.482 2.125 7.221 THROUGH 314 47.482 -2.125 7.221 
315 47.482 1.25 19.025 THROUGH 316 47.482 -1.25 19.025 
317 47.482 1.25 21.025 THROUGH 318 47.482 -1.25 21.025 
319 47.482 0.5 28.694 THROUGH 320 47.482 -0.5 28.482 
C 

C FWD SPAR, SHEAR WEB, (+45 TO +90) AND (-45 TO -90). 
C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 16,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 1 19 21 6 1 
GENERATE CONNECTS 4 22 24 6 1 
GENERATE CONNECTS 55 73 75 6 1 
GENERATE CONNECTS 58 76 78 6 1 
C 

C FWD SPAR, SHEAR WEB, +45 TO -45. 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 18,2 
QUADRILATERIAL PLATE TYPE 1 3 
GENERATE CONNECTS 19 31 33 6 1 
GENERATE CONNECTS 43 55 57 6 1 
GENERATE CONNECTS 22 34 36 6 1 
GENERATE CONNECTS 46 58 60 6 1 
GENERATE CONNECTS 31 43 44 6 1 
GENERATE CONNECTS 32 33 45 1 6 
GENERATE CONNECTS 35 34 46 1 6 
GENERATE CONNECTS 35 47 48 6 1 
LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 20,2 
QUADRILATERIAL PLATE TYPE 1 3 
GENERATE CONNECTS 79 86 89 6 2 
GENERATE CONNECTS 80 86 90 6 2 
CONNECT 39 81 79 33 
CONNECT 40 82 80 34 
CONNECT 39 81 83 45 
CONNECT 40 82 84 46 
C 

C FWD SPAR, SPAR CAP, +45 TO -45. 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 19,2 
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QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 45 57 58 6 1 
GENERATE CONNECTS 21 33 34 6 1 
LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 11,2 
QUADRILATERAL PLATE TYPE 1 3 
CONNECT 45 83 84 46 
CONNECT 33 79 80 34 
CONNECT 79 85 86 80 
CONNECT 83 89 90 84 
C 

C FWD SPAR, SPAR CAP, (+45 TO +90) AND (-45 TO -90). 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 19,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 57 75 76 6 1 
GENERATE CONNECTS 3 21 22 6 1 
C 

C FWD SPAR, INNER SPAR CAP, +90 TO -90. 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 12,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 2 74 77 6 3 
C 

C REAR SPAR, SHEAR WEB, (+45 TO +90) AND (-45 TO -90). 
C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 18,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 102 114 1 15 4 1 
GENERATE CONNECTS 101 113 116 4 3 
GENERATE CONNECTS 138 150 151 4 1 
GENERATE CONNECTS 137 149 152 4 3 
C 

C REAR SPAR, SHEAR WEB, +45 TO -45. 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 18,2 
QUADRILATERIAL PLATE TYPE 1 3 
GENERATE CONNECTS 130 138 139 4 1 
GENERATE CONNECTS 114 122 123 4 1 
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GENERATE CONNECTS 129 137 140 4 3 
GENERATE CONNECTS 1 13 121 124 4 3 
GENERATE CONNECTS 122 123 131 1 4 
GENERATE CONNECTS 121 124 132 3 4 
GENERATE CONNECTS 153 159 163 6 2 
GENERATE CONNECTS 154 160 164 6 2 
CONNECT 128 156 154 124 
CONNECT 127 155 157 131 
CONNECT 127 155 153 123 
CONNECT 128 156 158 132 
C 

C REAR SPAR, SPAR CAP, +45 TO -45. 

C 

LAMINATE DATABASE LAM. ARC 
MATERIAL COMPOSITE 15,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 115 123 124 4 1 
GENERATE CONNECTS 131 139 140 4 1 
GENERATE CONNECTS 113 137 138 4 1 
CONNECT 153 159 160 154 
CONNECT 157 163 164 158 
CONNECT 123 153 154 124 
CONNECT 131 157 158 132 
C 

C REAR SPAR, SPAR CAP, (+45 TO +90) AND (-45 TO -90). 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 19,2 
QUADRILATERAL PLATE TYPE 1 3 
GENERATE CONNECTS 101 113 1 14 4 1 
GENERATE CONNECTS 103 1 15 1 16 4 1 
GENERATE CONNECTS 137 149 150 4 1 
GENERATE CONNECTS 135 151 152 4 1 
C 

C RIGID WING STRUCTURE DESIGNED TO EXHIBIT MINIMAL 
C DEFORMATION WITHIN THE COMPONENT. 

C 

LAMINATE DATABASE LAM.ARC 
MATERIAL COMPOSITE 1,2 
QUADRILATERAL PLATE TYPE 1 3 
CONNECT 89 311 313 90 
CONNECT 85 312 314 86 
CONNECT 89 311 312 85 
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CONNECT 90 313 314 86 
CONNECT 163 317 318 159 
CONNECT 164 315 316 160 
CONNECT 163 164 315 317 
CONNECT 159 160 316 318 
CONNECT 163 307 319 317 
CONNECT 159 308 320 318 
CONNECT 90 164 315 313 
CONNECT 86 160316 314 
GENERATE CONNECT 307 308 320 1 12 
GENERATE CONNECT 301 302 310 1 8 
CONNECT 301 89 311 309 
CONNECT 302 85 312 310 
CONNECT 301 302 85 89 
CONNECT 85 86 90 89 
CONNECT 86 90 163 159 
CONNECT 163 159 160 164 
CONNECT 160 164 307 308 
GENERATE CONNECT 309 310 320 1 2 
END 



B. MSC/PAL2 LOAD FILE 

DISPLACEMENTS APPLIED 21 
ALL 0 1 73 74 72 1 
ALL 0 5 77 78 72 1 
ALLO 101 149 150 48 1 
DISPLACEMENTS APPLIED 1 
ALL 0 103 104 151 152 
ALL 0 3 4 75 76 

FORCES AND MOMENTS APPLIED 0 

FY 314 480 

SOLVE 

QUIT 
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C. MATLAB® PLY ANGLE OPTIMIZATION PROGRAM 

% ABD MATRIX OPTIMIZATION 
% 1.) Enter physical properties EEll through G12. 

% 2.) Select the desired QBxx values to be summed 
% by QBTOT. 

% 3.) Enter the correct ylabel value. 

EE11=1.97E7; 

EE22=1.2E6; 

U12=2.549E-1; 

U21=1.558E-2; 

G12=9.368E5; 

Q11=EE11/(1-U12*U21); 

Q22=EE22/( 1 -U 1 2*U21 ); 

Q12=U21*EE11/(1-U12*U21); 

Q66=G12; 
b=linspace(0,90,9 1 ); 

QBTOT=eye(l,91); 
for Z=l:l;91, 

B=b(Z)*3. 1415927/1 80; 

QBl 1=Q1 1 *(cos(B))M+2*(Q12+2*Q66)* 

(sin(B)^2)*(cos(B))/^2+Q22*(sin(B))M; 

QB22=Q1 1 *(sin(B))M+2*(Ql 2+2*Q66)* 

(sin(B)^2)*(cos(B))^2+Q22*(cos(B))M; 

QB 12=(Q1 l+Q22-4*Q66)*(sin(B))^2*(cos(B))^2+ 
Q12*((sin(B))M+(cos(B))M); 

QB66=(Q 1 1 +Q22-2*Q1 2-2*Q66)*(sin(B))^2* 
(cos(B))^2+Q66*((sin(B))M+(cos(B))M); 

QB 1 6=(Q1 1 -Q22-2*Q66)*sin(B)*(cos(B))^3+ 
(Q12-Q22+2*Q66)*(sin(B))^3*cos(B); 
QB26=(Qll-Q22-2*Q66)*(sin(B))^3*cos(B)+ 
(Q12-Q22+2*Q66)*sin(B)*(cos(B))^3; 

QBTOT(Z)=QB 1 1 +QB 1 2; 
end 



plot(b, QBTOT); xlabel(’Ply angle, degrees') ly label('Al 1+A12’):grid 
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APPENDIX C 



A. PLY PROPERTIES 





MATERIAL 


MATERIAL 


MATERIAL 




1 


2 


3 


El 


3.1 OOE+05 


7.568E+06 


1 .945E + 07 


E2 


2.067E+05 


1.994E+06 


1.193E+06 


G12 


2.314E+05 


8.433E+05 


9.135E + 05 


U12 


0.270E+00 


0.259E+00 


0.256E+00 


ALPHl 


2.220E-05 


2.01 9E-06 


8.564E-08 


ALPH2 


1.056E-05 


1 .288E-05 


1 .785E-05 


ALPH12 


O.OOOE+00 


O.OOOE+00 


O.OOOE + 00 


RHO 


6.714E-02 


7.1 OOE-02 


5.626E-02 


FTl 


O.OOOE+00 


2.980E+05 


2.500E + 05 


FT2 


O.OOOE+00 


1.500E+04 


1 .500E + 04 


FCl 


O.OOOE+00 


-1 .180E + 05 


-1 .600E + 05 


FC2 


O.OOOE+00 


-3.000E+04 


-3.000E + 04 


F12 


O.OOOE+00 


9.000E+03 


1.400E + 04 


FAW 


1.085E-04 


3.300E-02 


6.600E-02 


RC 


3.000E-01 


2.648E-01 


3.175E-01 


W 


7.143E-02 


2.000E-02 


2.000E-02 


TPLY 


2.309E-03 


6.322E-01 


1 .71 9E+00 


G13 


O.OOOE+00 


O.OOOE+00 


O.OOOE + 00 


G23 


O.OOOE+00 


O.OOOE+00 


O.OOOE + 00 


F13 


O.OOOE+00 


O.OOOE+00 


O.OOOE + 00 
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B. LAMINATE LAYUP 



1. Laminate layup 

LYR MAT THICK 

1 1 0.00225 

2 1 0.00225 

3 1 0.00225 

4 1 0.00225 

5 1 0.00225 

6 1 0.00225 

7 1 0.00225 

8 1 0.00225 

9 1 0.00225 

10 1 0.00225 

1 1 1 0.00225 

12 1 0.00225 



2. Laminate layup 

LYR MAT THICK 

1 1 0.0022 

2 1 0.0022 

3 1 0.0022 

4 1 0.0022 

5 1 0.0022 

6 1 0.0022 

7 1 0.0022 

8 1 0.0022 

9 1 0.0022 

10 1 0.0022 

1 1 1 0.0022 

12 1 0.0022 

13 1 0.0022 

14 1 0.0022 

15 1 0.0022 

16 1 0.0022 



location A 



ANGLE 


ID 


45.0 


12 


-45.0 


11 


45.0 


10 


-45.0 


9 


45.0 


8 


-45.0 


7 


-45.0 


6 


45.0 


5 


-45.0 


4 


45.0 


3 


-45.0 


2 


45.0 


1 



location 'B' and 'D' 



ANGLE 


ID 


45.0 


16 


-45.0 


15 


45.0 


14 


-45.0 


13 


45.0 


12 


-45.0 


11 


45.0 


10 


-45.0 


9 


-45.0 


8 


45.0 


7 


-45.0 


6 


45.0 


5 


-45.0 


4 


45.0 


3 


-45.0 


2 


45.0 


1 
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3. 


Laminate layup 


location 'C 




LYR 


MAT 


THICK 


ANGLE 


ID 


1 


1 


0.00225 


45.0 


20 


2 


1 


0.00225 


-45.0 


19 


3 


1 


0.00225 


45.0 


18 


4 


1 


0.00225 


-45.0 


17 


5 


1 


0.00225 


45.0 


16 


6 


1 


0.00225 


-45.0 


15 


7 


1 


0.00225 


45.0 


14 


8 


1 


0.00225 


-45.0 


13 


9 


1 


0.00225 


45.0 


12 


10 


1 


0.00225 


-45.0 


11 


11 


1 


0.00225 


-45.0 


10 


12 


1 


0.00225 


45.0 


9 


13 


1 


0.00225 


-45.0 


8 


14 


1 


0.00225 


45.0 


7 


15 


1 


0.00225 


-45.0 


6 


16 


1 


0.00225 


45.0 


5 


17 


1 


0.00225 


-45.0 


4 


18 


1 


0.00225 


45.0 


3 


19 


1 


0.00225 


-45.0 


2 


20 


1 


0.00225 


45.0 


1 



4. 


Laminate 


layup location 'AA' 


AND EE 


LYR 


MAT 


THICK 


ANGLE 


ID 


1 


2 


0.0070 


0.0 


6 


2 


2 


0.0070 


0.0 


5 


3 


2 


0.0070 


0.0 


4 


4 


2 


0.0070 


0.0 


3 


5 


3 


0.0120 


0.0 


2 


6 


3 


0.0120 


0.0 


1 
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5. 


Laminate layup 


location 


'BB' 


LYR 


MAT 


THICK 


ANGLE 


ID 


1 


2 


0.0070 


0.0 


11 


2 


2 


0.0070 


0.0 


10 


3 


2 


0.0070 


0.0 


9 


4 


2 


0.0070 


0.0 


8 


5 


2 


0.0070 


0.0 


7 


6 


2 


0.0070 


0.0 


6 


7 


2 


0.0070 


0.0 


5 


8 


2 


0.0070 


0.0 


4 


9 


3 


0.0120 


0.0 


3 


10 


3 


0.0120 


0.0 


2 


11 


3 


0.0120 


0.0 


1 



6. 


Laminate layup 


location 


•cc 


LYR 


MAT THICK 


ANGLE 


DD 


1 


2 0.0070 


0.0 


3 


2 


2 0.0070 


0.0 


2 


3 


2 0.0070 


0.0 


1 



7. 


Laminate layup 


location *DD* 


LYR 


MAT THICK 


ANGLE ED 


1 


2 0.0070 


0.0 4 


2 


2 0.0070 


0.0 3 


3 


2 0.0070 


0.0 2 


4 


3 0.0120 


0.0 1 
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LAMINATE PROPERTIES 
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1.427E+04 -2.385E+03 -U53E-06 -5.148E-05 9.711E-06 -1.459E-06 

-2.385E+03 1.427E+04 1.410E-06 9.711E-06 -5.148E-05 -1.359E-06 

-1.153E-06 1.410E-06 3.832E+03 -1.459E-06 -1.359E-06 -1.117E-05 
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Laminate 'AA* and 'EE* properties 



o o 
o o o 
-h + -f- 
PQ W W 
O O ON 

o o 
o o rN 

O O (N 



o o 
o o 
+ + 

W W PQ 
o o a^ 

O O (N 

o o o 

d d 



o 

+ 



m 

o 



o o 
o o 
+ + 
WWW 

O O ON 

o o m 

p p p 
d d in 



o o 
o o 

+ “h 

w w 
o o 
o o 
o o 

d d 



csi 

o 

I 

W 

cs 

r-- 

ON 



O 

W 

oo 

cn 

p 

d 



0 

X 

01 



w 



^ ^ o 
o o § 

+ + + 
W W ^ 
t =5 

^ ON § 

n ^ o 
r- <N 



o o 
o o 
+ + 
W W 
o o 
o o 
o o 

d d 



o o 
+ + 
w w 

r- oo 
O in 
CNJ in 

(N od 



o 

+ 

w 

in 

p 



o 

o 

+ 

w 

o 

o 

o 



in ^ O 
o o o 
+ + + 
WWW 
^ o 
^ o o 
oo (N O 

d (N d 



0 — 0 
O O O 
+ + + 
WWW 
^ — o 
ON o 
oo oo o 

d — ‘ d 



o o 
o o 
+ + 
W W 
o o 
o o 
o o 

d d 



00 



o o 
o o 
+ + 
W W 
o o 
o o 
o o 

d d 



^ m ^ § 

^ p p ^ 

5 w w w 

in 
cn 00 
— r- 



W 

I 

w 



o 
o 
o 
w d 



cn 




o 


(N 


o o 


in 


O 


O 


o 


O 


o o 


o 


+ 


+ 


+ 




+ + 


1 


W 


w 


W 


w 


w w 


w 


<N 


NO 


o 


in 


^ o 




r- 




o 


ON 


r- o 


NO 


p 


p 


o 


p 


oo o 


p 


d 


d 


d 


W 


d d 


d 



r-- 



o 

o 

o 

d 



o o 
o o 
+ + 
w w 
o o 
o o 
o o 

d d c4 



m 

o 

w 

ON 

ON 



r-- in g 
p p ^ 
WWW 
tn ON o 
NO cn o 

^ (N p 

d d 



S 

99? 

W w tu 
o >0 o 
m \o o 

O 

— ‘ 9 o 



2 9 o o 

^ S w w 

W 2 NO o 

Q 8 

T d d 



cn 

p 

W 

o 

(N 



cn S 

9 + 
w w 
o 
— o 

ON o 

w d 



o o 
o o 
+ + 
w w 
o o 
o o 
o o 

d d 



a:: 

Cu 

I 

00 



o o 
o o 
+ + 
w w 
o o 
o o 
o o 

d d 



o 

o 

+ 

W 

o 

o 

o 



o 

o 

+ 

w 

o 

o 

o 



Is 

UJ UJ 



99 

On q 

>9 »o 



tu H 

u z 

a: uJ 

o s 

tL, 3 

o ^ 

oo 

E- 

Z 

UJ 

u 

E 

UJ 

O 

U 

J 

< 

QC 

m 

X 

E- 



o 

o 

+ 

w 

o 

o 

o 

d 

II 

X 

< 

X 

w 

J 

< 



in 

O 

w 

o 

in 



< 

X 

w 

J 

< 



o 

o 

+ 

w 

o 

o 

o 

d 

II 

X 

S 

X 

w 

J 

< 

o 

w 

in 

r- 

cn 



X 

w 

w 

< 



m 

O 



NO 

o 



“s 

o ^ 

O P 

II 9 
X A 
< X 

ac 2 

£ X 

j ft. 



OO 

m 

(N 

d 

II 

X 

w 



in 

O 

+ 

w 

r-- 

p 

oo 

II 

X 

a 



p 

o ^ 
+ ^ 
W 

c<i 9 
m o 

II 

— ' Q/ 

sli < 

> CO 

W N 



+ o 
tu <N 

(SI ^ 
^ O 



X 

m 



3 G 



57 



Laminate 'BB' properties 
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Laminate 'CC properties 
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Laminate 'DD' properties 
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EX=1.019E+07 EY= 1.641 E+06 GXY= 8.680E+05 UXY=0.258 RH0= 2.166E-03 

TTH= 0.03300 ZBAR= 0.01650 



APPENDIX D 



A. TEST STAND CALCULATIONS 

Static analysis was performed on the components listed in Table Dl. 
Table Dl. COMPONENT ANALYSIS. 



Component 


Analysis 


Bearings 


Static load. 


Shaft 


Bending, shear, torsion. 


Extension 


Bending, shear, torsion. 


Bracket 


Bending, shear torsion. 



The force on each load-carrying member was calculated from the weight of the 
components and the forces applied during engine runs. A vector representing 
the weight of each component was place at its center of gravity as shown in 
Figure Dl. The engine loads were calculated for a 29 Hp motor producing 150 
lbs of static thrust. 

A strengths of materials approach was performed using the equations Dl 
through D8 [Ref. 15]. 



static load; 


If 


= 0 


(eq. Dl) 




Im 


= 0 


(eq. D2) 


shear: 


T 

max 


V 

= k — , solid beam 
A 


(eq. D3) 




X 

max 


= , at neutral axis 

It 


(eq. D4) 


bending: 


a 

max 


Me 

I 


(eq. D5) 
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torque: 



T c 

'^max = , circular shaft (eq. D6) 



2ta 



normal force: 



P 

a = — 
A 



, square beam (eq. D7) 



(eq. D8) 



The maximum stress components were calculated for each member and applied 
to Von Mises’ yield criterion (eq. D9) to compare to the material yield values. 



(c^xx - + (fjyy " )^ + (^^xx ~ + 6CJ„ + 



< — (eq. D9) 




Figure Dl. Test stand components. 
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A load applied to the engine produces an engine torque which creates torsion 
on the axial members of the test stand. Torque produce by the engine is a 
function of RPM and horsepower as related in equation DIO. The AROD uses a 
Herbranson engine rated at 29 Hp at 8000 RPM. 



-^J-(60=)*(l2i:) (eq.DlO) 

B. BEARINGS 

The bearings used were a Thompson PB-24 and a Thompson Super24 linear 
bearing. Description of these two bearings are listed is Table D2 [Ref. 15, p.l9]. 
The bearings are mounted in self-aligning pillow blocks 9 inches apart. The sum 
of the forces in the vertical direction was calculated using equation D1 to 
determine the static load on each bearing. 



Torque (in. •lbs. ) — 



Hp 

2 71 RPM 



Table D2. BEARING CHARACTERISTICS. 



Bore diameter, 
inches 


Bearing diameter, 
inches 


Length, 

inches 


Normal load rating, 
lbs. 


Maximum load 
rating, lbs. 


1.5 


2.375 


3.0 


1560 


2000 
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(eq. Dl) 



If 

F 



vertical 



= 0 



+ F +F +F +F +F =0 

Bonng A Bcanog B AROD Bneket Exuqcsoq SbiA 

Fbc^a +Fb«-^b - 49 lbs. - 12 lbs. - 11. 5 lbs. - 12 lbs. = 0 

Fb«^a +Fb^b = 84.5 lbs. 



Sm 




Bearing B 



= 0 



(eq. D2) 



A Beanng A AROD AROD Bncka Bnckei Excenaao Exteotion 



+ - 0 

,(9in.) + 49lbs.(45in.) + 12 lbs. (32 in.) + 11.5 lbs. (21 in.) 

+ 12 lbs. (Sin.) = 0 



Fb^sa = 314.5 lbs. (eq.Dll) 

The value obtained from the sum of the moments was substituted into the 
equation representing the sum of the forces: 

Fb^a +Fb.™sb = 84.5 lbs. 

Fb.^a = 314.5 lbs. (eq.Dll) 

F,.^, = 230 lbs. (eq. D12) 



The maximum static load from Table D2 was compared with the calculated 
static loads, (eq. Dl 1) and (eq. D12), for each bearing. As shown, both bearings 
provide an ample load margin under the maximum rated value. 

C. SHAFT 

The shaft is constmeted from hardened, highly polished steel described in 
Table D3. All forces generated on the test stand act on the shaft. Two groups of 
forces are present while the AROD is mounted on the test stand. When tests are 
not being performed, the weight of the AROD, bracket and extension apply a 
shear and bending load on the shaft. During engine tests an additional normal 
force due to thrust and torque due to horsepower are applied. 
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Table D3. SHAFT MATERIAL PI 


^lOPERTIES. 


Material 


Treatment 


Ou ,psi 


Gy ,psi 


Rockwell HN 


Brinell HN 


AISI 440C 


Case hardened 


285,000 


275,000 


C60 


610 



The maximum shear load in the shaft was determined by summing the 
vertical forces acting to the left of the shaft shown in Figure Dl. Denoting V as 
the shear force (eq. Dl) becomes: 

y F =0 

vertical ^ 



(eq. Dl) 



~ +F +F +V =0 

AJtOD Bncket Extcosioo Sb&A 

- 49ibs - 12 lbs. - ll.Sibs. + V 



Shaft 



0 



= 72.5 lbs. 

The weight of the shaft was neglected since its weight is supported by the 
bearings and not the shaft itself. 

The maximum stress due to the vertical shear component was calculated 
using equation D3, where k = 4/3 for a circular shaft. 

V 



= k 



A 

( 



(eq. D3) 



72.5 lbs. 
jr(0.75m.)' 



''max = 54.7 psi. 

A bending moment applied to the shaft results from the vertical forces 
applied to a moment arm originating at the left of bearing A. For these forces 
equation D2 becomes: 



Sm 



Shaft 



= 0 



(eq. D2 



^arod'^arod 0 

49(34.5) + 12(21.5) + 11.5(10.5) + = 0 

M,,. = 2069.25 in • lbs. 
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Since the shaft is symmetric the compressive and tensile stresses are of equal 
magnitude and opposite sign. The maximum value is computed using equation 
D5. 

o™. = ^ (eq. D5) 

2069.25 in. .lbs. * 0.75 m. 

0.2485 m.‘ 

= 6245.1 psi 

where I = ittr'' and c is the radius of the shaft. 

4 

The values for the engine were substituted into equation DIO. 

^ ^ 29hp* 550^^*60^*12.^ 

2 71 * SOOOrpm 
Tsh.ft = 228.5 in.*lbs. 

The shear stress due resulting from torsion was calculated using equation D6. 

t... ' ^ (eq. D6) 

(228.5 in..lbs.)(0.75 m.) 

(0.4970 m.") 

= 344.8 psi. 

where J is the polar moment of inertia, J = Itic" and c is the radius. 

The final force applied to the shaft is the normal force due to the static thrust 
of the ducted fan. For this calculation the maximum thrust output was estimated 
at 150 lbs. Substituting this value for P into equation D8 provides 



a = — (eq. D8) 

A 

150 lbs. 

71 (0.75 m.)‘ 
a = 84.9 psi. 
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The results of these maximum values of shear, bending, torque and normal 
force are depicted in Figure D2. 

Torque and normal force are constant throughout the shaft; however, 
bending and shear vary around the circumference. Bending is maximum on the 
upper and lower surface of the shaft and zero on the middle surface, while 
vertical shear is zero on the upper and lower surface and maximum on the 
middle surface. Two elements taken from the surface show the effects of these 
combined stresses. The stress values were applied to Von Mises’ yield criterion, 
(eq. D9), to compute the maximum stress seen by the shaft. 




Figure D2. Stresses on the shaft. 
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(eq. D9) 



- <^yyf + (Oyy ~ f + (c^xx ~ )' + 6CJy, + 6G„ + 60 , 

^(c^L+CyL+6axy)< Y 

V^L+3< <cj^. 

(a.) V84.9V. + 3(399.5 p.i)' < 275,000 ps. and (b.) V6330Vi + 3(344.8psif <275,000p« 

697,lps. < 275,000p« 6358.1psi < 275,000psi 



where Y is the yield strength from Table D2. Using Von Mises' yield criteria, 
the allowable load applied to the shaft produces stresses well below the yield 
point of the material. 

D. EXTENSION 

The extension was designed to locate the duct intake far enough away from 
the base of the stand to prevent any blockage of airflow. Since the extension was 
configured in a cantilever arrangement, it must support its own weight, that of 
the bracket and AROD, and any forces produced by the engine. The material 
selected for the extension was a low carbon steel, square tube. Its properties are 
listed in Table D4. 

Table D4. EXTENSION MATERIAL PROPERTIES. 



Material 


Dimensions 


Wall thickness 


Gtu .psi 


Gty xpsi 


CTcu Tsi 


CJsYxpsi 


AISI 1025 


2 in. X 2 in. 


0.125 in. 


55,000 


36,000 


36,000 


35,000 



The static analysis performed on the extension followed the same procedure 
as the analysis used for the shaft. The technique for calculating vertical shear 
and shear due to torsion varied slightly because of the different cross sectional 
area of the extension. 
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The vertical shear stress was calculated as 



^ ^vertical ~ ^ 

F + F + F 

^ AROD ^ Bncka ^ E: 




V 



EX(£DS10D 



72.5 lbs. 



(eq. Dl) 



X 



max 



T 



max 



YQ 

It 

(72.5 ibs.)(Q. 330 
(0.55l8m.‘)(0.25 in.) 
173.5 psi. 



(eq. D4) 



where Q= Ay is the first moment of area and 1= is the moment of 

inertia. 

The maximum bending moment was calculated at the right end of the 
extension where the longest moment arm occurs. The maximum normal stress 
became 



y =0 (eq. D2) 

Exiensioo x T. ^ 

AROd'^AROD Bncka^Bnckei Exteosioo ^ Exteosioo '^Exicmioo 

49lbs.(31in.) + 12lbs.(18m.) + 11.5 lbs. (7 in.) + M, = 0 

= 1815.5 in.^lbs. 



a 



max 



a 



max 



Me 

I 

(1815.5 m..lbs.)(lm.) 

(0.5518m.^) 
3290 psi. 



(eq. D5) 



where c is the distance from the neutral axis to the outer wall of the extension. 

The engine torque and the thrust were identical to the values derived for the 
shaft. The shear and normal stress were calculated using equation D7 and D8. 
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2 ta 



(eq. D7) 



(228. 5 in. •lbs.) 



2(0.125 in.)(3.5156u,^) 
260 psi. 



£ 

A 



(eq. D8) 



150 lbs. 



0.9375 in^ 

= 160 

where OC is the area bound by the centerline of the tube walls and A is the cross 
sectional area of the square tube. 

A combination of the stress values derived from maximum shear (a.) and 
maximum bending (b.) is shown in Figure D3. 



a.) 



T:xy - j 



Tc . 4V 
3 A 
= 260+ 173.5 
= 433.5 psi 

A I 
= 160 + 0 
= 160 psi 



b.) 



j 3 A 
= 260 + 0 
= 260 psi 



_z. Me 
‘Jxx - A + I 




=160 + 3290 
= 3450 psi 



Figure D3. Stresses on the extension. 
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These values also were substituted into Von Mises' yield criterion. 



^[(<^xx - <^yyf + (<^yy + («^xx " yz + + 6^ 

V<+3(J^y <^ys 



(eq. D9) 



(a.) Vl92%s. + 3(433.5 ps.)^ < 36,000 psi and (b.) V3450Vi + 3(260 psi)^ < 36,000p« 

775 psi < 36,000 psi 3479. 3 psi < 36, 000 psi 



E. BRACKET 

Forces from the AROD are supported by the four arms of the bracket. The 
bracket was configured in an "X" to allow the fuel tank vents to be located at the 
highest possible elevation. Figure D4 shows a front and side view of the bracket 
configuration. 




The material used for the bracket was aluminum square tubing. The 
material properties are given in Table D5. 
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Table D5. BRACKET MATERIAL PROPERTIES. 



Material 


Dimensions 


Wall thickness 


Gtu >psi 


Gty sPsi 


Gcu »psi 


Gsy ,psi 


6061-T6 


2 in. X 2 in. 


0.125 in. 


42,000 


37,000 


42,000 


20,000 



Each arm of the bracket reacts differently to the forces of thrust, weight, and 
engine torque. The forces are additive in some locations and cancel one another 
in others. A diagram of normal and shear stress was drawn to identify which 
locations would undergo the maximum stress condition. Figure D5 illustrates the 
stress combination on the compressive and tensile side of the bracket. 

Figure D5 identified the lower right arm of the bracket as experiencing the 
maximum combined stress condition. The values of normal and shear stress 
resulting from the forces on the bracket were calculated for the location 
immediately outboard from the extension flange were the greatest moment arm 
will occur. The resulting moment arm, 12 inches long, was used in the 
following calculations. 

1. Engine torque 

A rotational force is produced as the engine torque attempts to spin the 
bracket. This force produces a tangential shear in the rotational directional. A 
shearing stress (eq. D4) was obtained from the tangential shear. 



V 



engine torque 



Torque 

distance 



228. 5 in.*lbs. 



19 lbs. , total 



V 



torque 



X 



torque 



= 4.75 lbs. , per arm 
4 ^ 

VQ (4.75ibs.)(0.3301in.’) , , . 

— ^ r = 11.37psi. 

It (0.5518in.^)(0.25in.) 
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Figure D5. Bracket stress components. 



2. Engine thrust 

A bending moment on the bracket arms results from static thrust 
applying a normal force at the outer duct diameter. A normal stress (eq. D5) 
was calculated from this bending moment. 
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M 



total thxust 



— + ^arm “ 150 lbs,® 12 in — 1800 lb. •in. , tOtSJ 



M 



M 



thrust 



thrust 

= 450 in. lbs. , per arm 



Me (450 in . lbs.)(l in.) rco a 

— - = 652.4 psi. 

I 0.55 18 in.' 



3. Weight 

The weight of the AROD and bracket were reduced to a single force 
located 13.34 inches from the centerline between the tensile and compressive 
faces of the bracket. The resultant 61 -pound force was transformed into 
components perpendicular and parallel to the bracket arm. The component 
forces are shown as Fi and F2 in Figure D5. 

l^AROD l^Brackct 

— 49 lbs. + 12 lbs. 

Ftoul =61 lbs. 

F = Ii2!iL= 15.24 lbs. 

4 

F, = Fj = Fcos(45^) 

F, =Fj =10.78 lbs. 

The F2 component creates a shearing stress (eq. D7) by applying a 
twisting moment to the bracket arm. 

T = F2X 

= 10.78 lbs.* 13.34 m. 

T = 143.8 in-.lbs. 

T 

^ ~2ta 

1 43. 8 in. *lbs. 

~ 2(.125 in.)(3.5156in.^) 

tp^ =163.6 psi. 
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The Fi component creates a moment on the lower right bracket arm. 
This moment applies opposing stresses on the tensile and compressive faces of the 
bracket. The forces on each face was calculated through the summation of forces 
(eq. Dl) and summation of moments (eq. D2). 






vertical 



= 0 



F - F = -F 

lens. ^comp. 

=-10.78».. 



^ ^resultant “ ^ 

F„_(12.34..)-F,„,,(14.34„.) = 0 
F„= 1.1621 F,^,. 

Solving these equations simultaneously and substituting the forces into equation 
D8 provided the normal stress. The Gf 2 value was only calculated for the tensile 
face since the compressive normal stress will be less. 



F =66. 5 lbs. and F_. = 77.28 ibs. 

comp. tens. 




77.28 lbs. 

" 0.25 in. 

Op^ = 309.2 psi. 

The shear and normal stress from the torque, thrust and weight components 
were added and applied to Von Mises' yield criterion. 



— ^lorque — 11-37 psi + 163. 6 psi — 174.94 psi 

— Ojpjyjj + ~ 052.4 psi + 309.2 psi — 961.6 psi 



6 





\2 



+ 6o^, + 6a. 



-f 6a 



xy 




V^xx + 3< <CJ^, 



V96l.6%si + 3(174.94 psi)^ < 37,000 ps. 



1008.2 PS. <37,000 PS. 



(eq. D9) 
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